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ABSTRACT 


High resolution observations of 12 CO(2/'3), 13 CO(3. /, 9), and HCN (5."4) 

J— 1—0 in the galaxy merger Arp 299 (IC 694 and NGC 3690) show the line 
ratios vary dramatically across the system. The 12 CO/ 13 CO ratio is unusually 
large, 60 ± 15, at the IC 694 nucleus, where 12 CO emission is very strong, and 
much smaller, 10 ± 3, in the southern extended disk of that galaxy. Elsewhere, 
the 12 CO/ 13 CO line ratio is 5-20, typical of spiral galaxies. The line ratio 
variation in the overlap between the two galaxies is smaller, ranging from 10 ± 3 
in the east to 20 ± 4 in the west. 

The 12 CO/HCN line ratio also varies across Arp 299, although to a lesser 
degree. HCN emission is bright towards each galaxy nucleus and in the 
extranuclear region of active star formation; it was not detected in the IC 694 
disk, or the eastern part of the overlap region, leading to lower limits of 25 and 
20 respectively. By contrast, at the nuclei of IC 694 and NGC 3690 the ratios 
are 9 ± 1 and 14 ± 3 respectively. In the western part of the overlap region it is 
11 ±3. 

The large 12 CO/ 13 CO 1-0 intensity ratio at the nucleus of IC 694 can 
primarily be attributed to a low to moderate optical depth (t< 1) in the 12 CO 
1-0 line. These data support the hypothesis that unusually high 12 CO/ 13 CO 
line ratios (> 20) are associated with extremely compact molecular distributions 
in the nuclei of merging galaxies. Relative to 12 CO, 13 CO 1-0 is brightest in 
quiescent regions of low 12 CO surface brightness and weakest in starburst regions 
and the galactic nuclei. A medium consisting of dense (n = 10 4 — 10 5 cm -3 ) and 
warm ( T k > 50 K) gas will reproduce the extreme line ratios observed in the 
nucleus of IC 694, where the area filling factor must be at least 20%. 



Subject headings: galaxies: evolution — galaxies: individual (Arp 299) - 
galaxies: ISM — galaxies: starburst — radio lines: galaxies — radio lines: ISM 
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1. Introduction 

The inner kiloparsecs of starburst and interacting galaxies harbor stunning amounts 
of molecular gas, 10 9 — 1O 1O M 0 (e.g., Scoville et al. 1991; Bryant & Scoville 1996). In 
these environments, molecular clouds are subject to intense radiation fields, supernovae 
explosions, winds from newborn hot stars, strong tidal forces, and gas surface densities 
several order of magnitudes higher than in the Milky Way disk. These are also extremely 
active star formation sites. Knowledge of the physical conditions and structure of molecular 
gas in interacting systems is essential to understand the starburst activity and its role in 
galaxy evolution. 

Arp 299 is an IR-luminous (Ljr & 8 x 10 11 L 0 ) merger system of two galaxies, IC 694 
and NGC 3690. Strong 12 CO emission has been detected from the nuclei of IC 694 and 
NGC 3690 and from the interface between the two galaxies (Solomon & Sage 1988; Casoli 
et al. 1989; Sargent et al. 1987; Sargent & Scoville 1991). The two nuclei, as well as the 
western overlap region, currently harbor intense star formation activity (c.f., Gehrz et al. 
1983; Baan & Haschick 1990). Furthermore, the nucleus of IC 694 is a flat-spectrum radio 
source, and may be an AGN (Gehrz et al. 1983). 

Lower resolution (single dish) observations reveal an unusually large 12 CO/ 13 CO 

1- 0 line intensity ratio, >20 in Arp 299 (Aalto et al. 1991; Casoli et al. 1992). These 
observations left it unclear whether this is due to weak 13 CO in the whole system or to a 
varying 12 CO/ 13 CO 1-0 line ratio. Observations at 20" and 11" resolution by Casoli et al. 
(1992) suggest little variation in the 12 CO/ 13 CO 1-0 line ratio and none in the 12 CO/ 13 CO 

2- 1 line ratio. In contrast, Aalto et al. (1995) note substantial variations at 28" resolution 
in the 12 CO/ 13 CO 2-1 ratio: the ratio is about 17 in IC 694; close to 9 in NGC 3690; and 
about 7 in the interface region between the two disks. Solomon et al. (1992) detected 
bright HCN emission in 28" maps of Arp 299. They measured 12 CO/HCN ratios of 11 in 



IC 694 and 13 in the interface region. 


2. Observations and Results 

Aperture synthesis 12 CO, 13 CO, and HCN 1-0 mapping of Arp 299 was carried out 
with the Owens Valley Radio Observatory (OVRO) millimeter array between March, 
1995, and February, 1996. SIS receivers on the six 10.4 m telccopes provided typical 
system temperatures (SSB) of 600 K, 450 K, and 350 K for 12 CO, 13 CO, and HCN. 
Quasars 1150+497 and 0917+449 were used for phase calibration and Uranus and Neptune 
for absolute flux calibraton. The synthesised beams are 2 ."5 x 2 ."2 for 12 CO (uniform 
weighting), 4."3 x 3."6 for 13 CO (natural weighting), and 5."6 x 5."3 for HCN (natural 
weighting). At 2.6 mm wavelength with 2."3 resolution, a brightness temperature of 1 K 
corresponds to 57 mJy beam -1 . The digital correlator, centered at Ulsr = 3100 km s -1 , 
provided a total velocity coverage of 1123 km s -1 for 12 CO, 1175 km s -1 for 13 CO, and 
1407 km s -1 for HCN. Data were binned to 4 MHz resolution, corresponding to 10 km 
s -1 for 12 CO and 13 CO and 13 km s -1 for HCN. At 110 GHz an unresolved continuum 
source of 17 ± 2 mJy was detected in the center of IC 694. Continuum emission was also 
detected in NGC 3690 (5 + 2 mJy) and the overlap region (9 + 2 mJy). We subtracted this 
continuum emission from the line emission before maps were made. 

The main structures found by Sargent & Scoville (1991) with the three telescopes 
array are recovered in our new 12 CO map (Figure la), but our increased uv coverage enable 
improved deconvolution. Lower surface brightness emission (A2), possibly a molecular disk 
or bar, extends 10-15" southeast of the IC 694 nucleus (Al), coincident with the remnant 
optical disk. The center of NGC 3690 (Bl) is also surrounded by weak, extended emission 
(B2 and B3) with a somewhat S-shaped morphology. Where the two galaxies overlap, three 
distinctive clumps (Cl, C2, and C3) can now be discerned. Weaker, extended emission is 
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also recovered better in our new 12 CO map and clumpy structures can be distinguished at 
the center of the map (F). These structures appear to connect the major components A, B, 
and C. In addition, there is a clump (D) north of the main features with systemic velocity 
~ 3280 km s^ 1 , which is beyond the bandwidth of the earlier OVRO data. 

Bright 6 cm radio continuum peaks at the nuclei of IC 694 and NGC 3690 (Gehrz 
et al. 1983; Condon et al. 1991) coincide with the 12 CO peaks to within the estimated 
positional uncertainty (0."5). There is also reasonable positional agreement between the two 
brightest 12 CO clumps in the overlap region (C) and two additional 6 cm radio continuum 
peaks: Cl and the western radio continuum peak are also coincident within 0. w 5, although 
the discrepancy between C2 and the eastern radio continuum peak is somewhat larger, 

A a = 0."8. Weaker radio continuum emission at 18 cm and 6 cm is spatially coincident 
with the extended 12 CO emission in regions A2, B2, B3, and F (Baan & Haschick 1990; 
Gehrz et al. 1983). 

The velocity held image (Figure 2, Plate 1) reveals a monotonic shift from the 
blueshifted emission from A2 in IC 694 to redshifted emission from the overlap region C and 
from region D. Velocity gradients within C and D are small. The velocity held of NGC 3690 
is complicated by a double-peaked emission structure of B2, and both B2 and B3 appear 
blueshifted relative to the center, Bl. 

The nuclei of both IC 694 (Al) and NGC 3690 (Bl) remain unresolved by our 2."2 
synthesised beam. A two-dimensional Gaussian ht to the nuclear emission of IC 694 and 
NGC 3690 yields upper limits to the source diameters of l."4 and l."5, corresponding to 
radii of 140 and 150 pc, respectively (for D =42 Mpc). This implies a lower limit of 18 K 
to the 12 CO brightness temperature in the IC 694 nucleus, and therefore, the cloud filling 
factor is quite high in the inner 200 pc. Even if the intrinsic brightness temperature is as 
high as 100 K, the surface filling factor of clouds is still almost 20%. The three features in 
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the overlap region are resolved, with sizes of 792 x 322 pc (Cl), 876 x 428 pc (C2), and 
684 x 456 pc (C3). Derived properties for all designated regions are presented in Table 1. 

The total molecular mass of Arp 299, estimated from M(H 2 ) = 1.2 x 10 4 SAr;D 2 M 0 
(D is the distance in Mpc; SAv is the integrated 12 CO 1-0 line flux in Jy K kms -1 ), is 
7.5 x 10 9 M 0 , 87% of the value estimated by Solomon and Sage (1988) from their single 
dish map. This formula corresponds to 1V(H 2 )//( 12 C0)=3 x 10 2 ° cm -2 (K km s -1 ) _1 
(Sargent & Scovillc 1991), a standard Galactic 12 CO luminosity to H 2 mass ratio. The 
conversion factor may vary, however, across Arp 299, since the line ratios indicate different 
cloud properties. 

The most striking feature of the 13 CO 1-0 map (Figure lb) is the absence of strong 
emission at the nucleus of IC 69f (Al). Emission is detected in the A2 region of IC 694, 
in the overlap region C, and in NGC 3690. Unlike Casoli et al. (1992) we End significant 
variation in the 12 CO/ 13 CO 1-0 line ratio across Arp 299. ^From the very high value of 
60 ± 15 at the IC 694 nucleus (Al) the ratio drops to 10 ± 3 in the A2 region. The global 
ratio for NGC 3690 is 13 ± 2, but there is an indication that the ratio is higher in the 
nucleus, Bl, than in B2 or B3. The ratio for the east (C3) and west (Cl) overlap region is 
11 ± 3 and 20 ± 4 respectively. The spatial correlation between the 13 CO emission and radio 
continuum emission peaks is also poor. 

As for 12 CO, the HCN map (Figure lc) is dominated by a peak at the center of IC 694 
(Al), where the 12 CO/HCN line ratio is 9 ± 1. Emission is also detected at the nucleus 
of NGC 3690, where the ratio is 14 ± 3, and in the western (Cl) overlap region, where 
12 CO/HCN= 11 ± 3. There is very little HCN in regions A2 and C3. The HCN and 13 CO 
line emission peaks appear anticorrelated. 



3. The molecular line-ratios 


As a system, Arp 299 is not deficient in 13 CO. In all regions but the core of IC 694, we 
observe 12 CO/ 13 CO line ratios typical of star forming regions in other galaxies (e.g., Aalto 
et al. 1995, 1991; Young & Sanders 1986; Rickard & Blitz 1985). Relative to 12 CO, 13 CO 
1-0 is brightest in quiescent regions of low 12 CO surface brightness, and weak in starburst 
regions and galactic nuclei. In contrast, HCN 1-0, like 12 CO, is bright in the two galaxy 
centers and in regions of active star formation. These line ratio variations are, most likely, 
caused by differences in line excitation. Our results support the suggestion (Aalto et al. 
1995) that unusually high (i.e. > 20) 12 CO/ 13 CO 1-0 line ratios tend to be associated 
with extremely compact molecular distributions centered on the nuclei of merging galaxies 
and are primarily due to a small or moderate optical depth, r<l, in the 12 CO 1-0 line. 
High ambient pressures, strong tidal forces and ongoing starburst or AGN activity lead to 
substantial changes in cloud structure and physical conditions. Bryant (1996) also finds 
that HCN 1-0 is bright in regions of high 12 CO surface brightness in merging galaxies. 


3.1. IC 694 

The line ratio variation from the IC 694 disk, where 12 CO/ 13 CO = 10 and 12 CO/HCN 
t 25, to its nucleus, where the corresponding values are 60 and 9, reflects a dramatic change 
in cloud properties. 

The bright HCN line accompanied by relatively weak 13 CO emission (HCN/ 13 CO = 

7 ± 2) implies a population of unusually dense and warm clouds. The HCN 1-0 strength 
implies densities n> 10 4 cm -3 if the HCN excitation is dominated by collisions with H 2 . 

It is also likely that the density is <10 5 , so that most of the HCN population will remain 
in the lower levels and r 10 > 1. At these densities, the 12 CO and 13 CO 1-0 transitions 
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are thermalised. If the kinetic temperature is also high, the lower levels may become 
significantly depopulated, effectively reducing the optical depth of the 1-0 line. Then the 
13 CO 2 1 /10 line ratio should be > 1. Comparing the single dish 13 CO 2-1 flux (Aalto et 
al. 1995) with our 13 CO 1-0 flux from the nucleus of IC 694, we estimate 13 CO 2—1 /1—0 
>2, implying that the gas temperature is high, > 50 K. Although the single dish beam was 
large, 28", the bulk (>65%) of the 12 CO 2-1 emission within the beam originates in the 
nucleus of IC 694 (Al). 

Since the lower transitions of 12 CO and 13 CO appear thermalised, LTE can be assumed 
to estimate the 12 CO column density, ./V( 12 CO). We find it unlikely that the 12 CO/ 13 CO 
abundance ratio is extremely high (section 3.3). Therefore, a high 12 CO/ 13 CO 1-0 line 
ratio indicates a low to moderate optical depth (t<1) in the 12 CO 1-0 line. The high 
intrinsic 12 CO 1-0 brightness temperature makes r << 1 unlikely, and we therefore 
assume ti 0 ( 12 CO) ~ 1. The optical depth of the 1-0 transition can be expressed as 
r 10 ( 12 CO) » 3.9 x 10~ 15 AT( 12 CO)(1 - e~ 5 - 53 / T ™)/T ex AV. For a temperature T ex =100 K, 
line width AV = 5 km s _1 , and ri 0 =l, the 12 CO column density ./V( 12 CO) = 2 x 10 18 
cm -2 (per cloud) and the resulting brightness temperature Tg( 12 CO) ~ 60 K. For a 
density of n = 10 4 cm -3 and a 12 CO abundance, [ 12 CO/H 2 ]=5 x 10~ 5 , the cloud radius 
r = A^( 12 CO)/2a;( 12 CO)n(H2)=0.7 pc — not unlike cores or clumps within Giant Molecular 
Clouds in the Galaxy. At this gas density, HCN is not thermalised so T ex (HCN) for the 
1-0 line will be significantly lower than that for 12 CO. If the abundance ratio [ 12 CO/HCN] 
ps 10 3 , then r 10 (HCN) is an order of magnitude higher than r 10 ( 12 CO). The resulting 
Tb(HCN) is sufficiently high to account for a 12 CO/HCN line ratio of 9. 

Above, we infer a clumpy molecular medium because we chose a AV = 5 km s” 1 , 
yielding small clouds with r < 1 pc (cf Aalto et al. 1991). We can not, however, exclude 
significantly larger AV which would indicate a continous, non-cloudy structure — perhaps 
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even a smooth, rotating disk. A third alternative is a molecular ISM consisting of dense 
clumps surrounded by diffuse, non-cloudy molecular gas (e.g. Aalto et al. 1994,1995; 
Dahmen et al. 1996). 


3.2. The overlap region C and NGC 3690 

Molecular line ratio variations are also seen within the overlap region C, albeit smaller 
than those within IC 694. The weakest 13 CO and strongest HCN 1-0 emission, relative to 
12 CO, is found in Cl, the location of the brightest Ha emission in Arp 299 (Gehrz et al. 
1983). Continuum emission at 3.4 /irn and 10/jrn also peaks close to the Cl clump. Thus, 
Cl appears to be currently the most active starforming region within C. We suggest that 
the observed molecular line ratio gradients are the result of a temperature and/or density 
gradient across the overlap region. The 12 CO/ 13 CO 1-0 line ratio in Cl is considerably 
lower than in the nucleus of IC 694, perhaps because Cl is an extranuclear starburst. 

Unlike 12 CO, the 13 CO emission peak is not connected with the nucleus of NGC 3690 
(Bl). The radial excitation gradient is similar, therefore, to that of IC 694, with the 
highest 12 CO/ 13 CO line ratio in the central region. The nucleus of NGC 3690, with intense 
associated Ha emission, is a site of starburst activity (Gehrz et al. 1983). 


3.3. Molecular abundances 

It has been suggested that high 12 CO/ 13 CO intensity ratios in mergers are caused by 
unusually high isotopic abundance ratios in molecular clouds with optically thick 12 CO 
1-0 lines. An influx of very low metallicity gas from the outer disk of the galaxies is the 
proposed cause of such an extreme abundance ratio (e.g. Casoli et al. 1992; Henkel & 
Mauersberger 1993). Since the 12 CO/ 13 CO intensity ratio is normal in the outskirts of 
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IC 694 and NGC 3690 this scenario is unlikely to be the explanation for Arp 299. Instead, 
the measured line ratio variations most likely indicate differences in the line excitation and 
gas properties in different parts of the system. 

The observed 12 CO/ 13 CO 1-0 line ratio is, however, a lower limit to the 12 CO/ 13 CO 
abundance ratio in the emitting region and for Tio( 12 CO) <1 this implies an abundance 
ratio not much greater than 60 in the center of IC 694. This value is typical for GMCs 
in the Galactic disk, but higher than in the inner region of our Galaxy, where the ratio is 
« 25 (Langer & Penzias 1990). Perhaps the ISM in the nucleus of IC 694 recently arrived 
from the disk of the galaxy. In this case, the difference in line ratio between A2 and Al is 
solely caused by a dramatic change in mean optical depth of the 12 CO line. On the other 
hand, selective photodissociation of 13 CO by a starburst and/or an AGN may change the 
isotopic abundance ratio. A young nuclear starburst may also produce extra 12 C and thus 
temporarily increase the 12 C/ 13 C abundance ratio (e.g., Henkel & Mauersberger 1993). 


4. Conclusions 

The 12 CO/ 13 CO 1-0 line ratio varies dramatically within Arp 299, from 60 ± 15 at the 
nucleus of IC 694 to 5-10 in its disk and in the eastern and north interface regions (C3 and 
D). The 13 CO 1-0 brightness, relative to 12 CO, is high in quiescent regions of low 12 CO 
surface brightness, and low in starburst regions and galactic nuclei. In contrast, HCN 1-0 
is bright in the two galaxy centers and in the active extranuclear star formation region. 
The 12 CO/HCN 1-0 is 9 ± 1 at the nucleus of IC 694, 14 ± 3 for NGC 3690 and 11 ± 3 for 
the extranuclear starburst region Cl. Unusually high 12 CO/ 13 CO line ratios (> 20) appear 
to be associated with extremely compact molecular distributions in the nuclei of merging 
galaxies (cf. Aalto et al. 1995). 
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The large 12 CO/ 13 CO 1-0 intensity ratio at the nucleus of IC 694 can be attributed 
to low to moderate optical depth (t< 1) in the 12 CO 1-0 line, possibly combined with 
abundance effects. A medium consisting of dense (n = 10 4 — 10 5 cm -3 ), warm (7b > 50 K) 
gas is consistent with the observations. 


We thank Peter Bryant for helpful discussions and suggestions. The OVRO mm array is 
supported in part by NSF grant AST 9314079 and the K.T. and E.L. Norris Foundation. 
The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Fig. 1.— (a) 12 CO 1-0 moment map (2800-3400 kms -1 ). Contour levels are 0.9, 1.8, 3.6, 5.4, 
7.2, 9, 10.8, 12.6, 14.4, 16.2, and 18 Jy beam -1 km s -1 . Gray scale levels range from 5 to 
75 Jy beam -1 km s -1 . The peak flux is 97 Jy beam -1 km s -1 on the nucleus of IC 694. 
Crosses mark the 6 cm radio continuum positions (Gehrz et al. 1983). (b) The 13 CO 1-0 
moment map. Contour levels are 0.4, 0.8, 1.2, 1.6, 2.0 Jy beam -1 km s -1 . The peak flux is 
2.3 Jy beam -1 km s -1 on the disk of 1C 694.(c)HCN 1-0 moment map. Contour levels are 
0.6, 1.2, 1.8, 2.4, 3.0, 3.6, 4.2, 4.8, 5.4 Jy beam -1 km s -1 . The peak flux is 7.1 Jy beam -1 
km s -1 at the nucleus of IC 694. 

Fig. 2.— The 12 CO 1-0 velocity field. The grayscale range from 2800 (light) to 3300 (dark) 
kms -1 , the contours from 2800 to 3200 kms -1 . 
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Table 1. Arp 299: Source properties 


Region 

Size a 

[arcsec] 

12 CO l-i 

Peak 

Jy kms -1 

beam -1 

0 Intensity 

Integrated 

[Jy kms -1 ] 

AVfwhm 

[kms -1 ] 

12 CO 1-0 

12 CO 1-0 

13 CO 1-0 HCN 1-0 

(note c) 

IC694: nucleus (Al) 

<1.4 

97 

122 

350 

60 ±15 

9 ± 1 

IC 694: disk (A2) 

17 x 6 b 

12 

120 

65 

10 ±3 

> 25 (3c) 

NGC3690: nucleus (Bl) 

<1.5 

29 

36 

260 



NGC3690: total 



53 


13 ±2 

14 ±3 

Overlap east (C3) 

3.4 x 2.2 

9 

15 

60 

10 ±3 

> 20 (3a) 

Overlap west (C2) 

4.3 x 2.1 

20 

40 

80 



Overlap west (Cl) 

3.9 x 2.0 

21 

47 

60 

20 ±4 

11 ±3 

Region D 



4 

80 

5 ± 3 

> 10 

Region F 
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40-80 




a FWHM from two dimensional Gaussian fits using the AIPS task IMFIT. 
b An estimate of the size of the SE disk of IC 694 — no Gaussian fit possible. 

c All ratios are in terms of integrated brightness temperatures. Uncertainties include thermal noise 
only. The 12 CO map was smoothed to the resolution of the 13 CO or HCN map before line ratios were 
constructed. 
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